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Abstract 

This  paper  focuses  on  designing  thermoelectric  generator  (TEG)  via  numerical  method.  Thermoelectric  module  in 
this  study  was  made  from  n-type  CaMn03  and  p-type  Ca3Co409.  The  couple  field  equations  were  formulated  into  a 
system  of  algebraic  equations  by  means  of  finite  element  method  (FEM).  The  simulation  procedure  was  verified  by 
comparing  with  experimental  results,  which  showed  a  very  good  agreement  between  the  two  outcomes.  The 
parametric  design  exhibited  that  output  voltages  vary  linearly  with  height  of  TEG,  until  reaching  some  limiting  value, 
the  linear  relationship  was  terminated.  After  that  specific  height,  increasing  the  height  of  TEG  no  longer  gains  notable 
output  voltage. 
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1.  Introduction 

The  development  of  materials  science  with  good  thermoelectric  properties  made  achievable  the 
fabrication  of  thermoelectric  generators  (TEG)  [1-3].  Being  sophisticated  devices,  thermoelectric  systems 
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offer  some  unique  advantages,  such  as  high  reliability,  long  life,  small-size,  environmentally-green  and 
no-vibrations  generator,  and  can  be  used  in  a  wide  temperature  range.  Thermoelectric  generators  are 
based  on  the  Seebeck  effect  for  power  generation.  If  a  steady  temperature  gradient  is  applied  along  a 
conducting  sample,  the  initially  uniform  charge  carrier  distribution  is  disturbed  as  the  free  carriers  located 
at  the  high-temperature  end  diffuse  to  the  low-temperature  end.  This  results  in  the  generation  of  a  back 
electromotive  force  which  opposed  any  further  diffusion  current.  The  open-circuit  voltage  when  no 
current  flows  are  the  Seebeck  voltage.  When  the  junctions  of  a  circuit  formed  from  two  dissimilar 
conductors  (n-  and  p-type  semiconductors)  connected  electrically  in  series  but  thermally  in  parallel  are 
maintained  at  different  temperatures  rhot  =  Tx  and  TcooX  =  T2  the  open-circuit  voltage  V  is  developed.  The 
thermoelectric  generator  is  a  standardized  device  consisting  of  p-  and  n-type  legs  connected  electrically  in 
series  and  thermally  in  parallel,  and  bonded  to  a  ceramic  plate  on  each  side.  The  modules  are  fabricated  in 
a  great  variety  of  sizes,  shapes  and  number  of  thermoelectric  couples  and  can  operate  in  a  wide  range  of 
temperature  gradient.  When  a  temperature  gradient  is  applied  across  the  thermoelectric  generator,  the  heat 
absorbed  at  the  hot  junction  and  cold  side  will  generate  a  current  through  the  circuit  and  deliver  electrical 
power  to  the  load  resistance.  This  study  aims  to  find  out  the  effect  of  geometric  design  to  the 
thermoelectric  generator.  The  typical  models  of  square-cross-section  area  TEG  associated  with  various 
heights  are  under  investigation.  By  assigning  area  to  constant,  but  increase  height  value,  we  could  come 
upon  the  optimal  height  which  issue  highest  power  output.  This  process  is  very  helpful  to  design  real 
thermoelectric  devices  for  commercial  products. 


2.  Computational  Details 


The  governing  equations  to  describe  behaviour  of  thermoelectric  material  are  governed  by  the  couple 
equations  of  heat  transfer  and  continuity  of  current  density  phenomena  as  follow  [4]: 

pc^  +  V-q=Q  (1) 

ot 


V- 


dE_ 

dt 


+  V-J  =  0 


(2) 


Where  T  is  temperature,  E  is  electric  field,  q  is  heat  flux,  J  is  electrical  current  density,  Q  is  internal 


heat  generator,  t  is  time,  p  is  density,  c  is  heat  capacity  and  s  is  electric  permittivity.  The  current  density 
J  is  generated  by  a  coupling  of  reversible  Seebeck  effect  and  irreversible  Joule  effect  can  be  written  as 


J  -  <jE  -  crSV T  (3) 

The  heat  flux  q  is  generated  by  couple  effect  of  reversible  Peltier  and  irreversible  Fourier  effect  is 
defined  as 


q=I\J-iNT  (4) 

where  El  is  Peltier  coefficient,  k  is  thermal  conductivity,  a  is  electrical  conductivity  and  S  is  Seebeck 
coefficient. 

The  electric  field  E  can  be  derived  from  an  electric  scalar  potential  (j)  as 

E  =  -V(/)  (5) 

The  two  coupled  governing  equations  (1)  and  (2)  are  then  transformed  into  finite  element  equations  by 
approximating  the  primitive  physical  unknowns,  temperature  T  and  electrical  potential  (j),  into 
interpolation  functions  and  value  of  nodal  known  on  a  element  by 
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T  =  [N]{Te}  (6) 

< *  =  [*]{&}  (7) 

where  Te  is  vector  of  nodal  temperature,  </)Q  is  vector  of  nodal  electrical  potential  and  N  is  element  shape 
functions.  After  lengthy  manipulations  base  on  the  Galerkin  weighting  scheme,  the  differential  equations 
become  algebraic  finite  element  equations  [5-7], 


CT  0 
0  CE 


The  global  matrix  equation  is  assembled  from  the  individual  finite  element  equations.  The  solution 
yields  temperatures,  Te ,  and  electric  potentials,  <f>e,  at  unconstrained  nodes,  or  reactions  in  the  form  of  heat 
flow  rate  and  electric  current  at  nodes  with  imposed  temperature  and  electric  potential  respectively 

Thermoelectric  module  is  made  of  CaMn03  and  Ca3Co409  for  n-and  p-type  leg,  respectively.  The 
module  is  fabricated  by  connect  n-and  p-cells  with  thin  copper  plates  and  superimposed  alumina  plates  on 
top  and  bottom  sides  of  its  structure,  as  shown  in  Fig.  1(a)  and  1(b).  To  establish  the  electrical  contact 
between  the  cells,  the  copper  and  alumina  plate  are  served  as  substrate  and  supporting  floor.  The  n-and  p- 
leg  was  similar  dimension  of  3x3x4  mm3.  Notice  that  the  height  of  n-and  p-leg  are  4  mm.  The  copper 
plate  and  alumina  plate  have  size  of  5x8x0. 3  mm3  and  5x10x0.7  mm3,  respectively. 
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Fig.  1.  (a)  &  (b)  Fabrication  of  a  single  module  TEG,  (c)  &  (d)  The  laboratory  equipment  set-up  for  measure  output  voltage  obtained 
from  differential  temperature  on  the  TEG 
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Fig.  2.  Geometric  detail  of  TEG 
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Table  1.  Material  properties 


Materials 

S  [pV/K] 

□  [□•m] 

|  [W/m-K] 

Ca3Co409 

147.31 

9.78xl0*5 

8.55 

CaMn03 

-480.23 

12.50xl0'5 

9.79 

Copper 

6.50 

1.69x10  s 

350 

Alumina 

- 

l.OOxlO14 

40 

Thermoelectric  structure  was  activated  by  differential  temperature  varying  in  range  0-40  K.  The  output 
voltage  generated  was  measured  by  the  instruments  as  presented  in  Fig.  lc  and  Id.  As  visualized  by  the 
plotting  results  in  Fig.  4(a),  output  voltage  related  to  driven  temperature  as  linear  function  exhibit  by  a 
straight  line. 

3.  Results  and  Discussion 

Numerical  simulation  is  carried  out  by  using  a  finite  element  package  ANSYS.  The  procedure  for 
doing  a  thermoelectricity  analysis  consists  of  following  main  steps:  create  the  physics  environment,  build 
and  mesh  the  model  and  assign  physics  attributes  to  each  region  within  the  model,  apply  boundary 
conditions  and  loads  excitation,  obtain  the  solution,  review  the  results. 


(a)  (b)  (c) 


Fig.  3.  (a)  Computational  mesh,  (b)  Temperature  distribution,  (c)  Voltage  distribution 


Simple  thermoelectric  generator  geometry  consists  of  n-type  CaMn03  and  p-type  Ca3Co409. 
Geometric  detail  of  TEG  is  shown  in  Fig.  2.  Number  encompassed  by  a  little  circle  prescribed  material 
type:  1  =  n-type  leg,  2  =  p-type  leg,  3  =  copper,  4  =  alumina.  The  material  properties  for  the  calculations 
with  temperature  independent  values  are  shown  in  Table  1.  In  the  present  application  for  modelling  the 
electric  and  thermal  fields  the  second  order  brick  element  was  chosen.  There  are  20  nodes  per  element. 
Each  node  consists  of  2  unknowns,  temperature  and  voltage. 
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(a)  (b) 

Fig.  4.  A  comparison  of  output  voltage  between  the  experiment  and  FEM  simulation  and  output  voltage  vs.  height  of  TEG 


First  of  all,  we  verify  the  computational  procedure  by  comparing  numerical  results  with  the 
experiment  in  section  III.  Fig.  4(a)  shows  good  agreement  between  results  from  laboratory  test  and  FEM 
solution.  Variation  of  the  voltage  output  varies  against  differential  temperature  related  in  linear 
relationship.  The  more  differential  temperature  results  in  the  more  voltage  output.  The  maximum  voltage 
23.76  mV  was  obtained  when  applied  40  K  differential  temperatures.  Distribution  of  temperature  and 
voltage  also  are  show  in  Fig.  3(b)  and  3(c),  respectively. 

To  study  effect  of  geometric  design,  we  characterize  thermoelectric  generator  into  various  height.  The 
height  of  the  module  varies  between  0-10  cm.  The  results  of  output  vary  against  height  were  plotted  in 
Fig.  4(b).  From  the  figure  we  obviously  found  that  there  are  3  regions  of  output  voltage.  Firstly,  when 
height  varies  increasingly  from  0  to  6  mm,  output  voltage  has  linear  relationship  with  height.  Secondly, 
when  height  varies  increasingly  from  6  mm  to  40  mm,  output  voltage  still  increasing  but  the  slope  of  the 
curve  was  drop.  Finally,  when  the  height  varies  is  increasing  from  40  mm  to  100  mm,  output  voltage  do 
not  increasing  but  exhibit  saturation  state.  The  numerical  results  suggested  that  we  can  use  Fig.  4(b)  in  a 
design  process  of  TEG.  This  guideline  leads  to  optimum  design.  The  amount  of  material  use  in 
fabrication  of  TEG  could  be  economical  by  using  prior  results  from  computational  simulation. 

4.  Conclusions 

In  this  paper,  we  use  numerical  technique  to  study  behaviour  of  thermoelectric  generator.  The 
governing  equations  were  transformed  into  algebraic  equations  via  FEM  discretization  technique.  In  this 
study,  we  analyzed  a  single  module  of  thermoelectric  generator  which  made  from  n-type  CaMn03  and  p- 
type  Ca3Co409.  Laboratory  experiment  was  performed  by  using  TEG  with  4  mm.  Height.  The  results 
from  experimental  test  shown  that  output  voltage  vary  linearly  with  differential  temperature.  Simulation 
results  from  FEM  were  closely  agreed  with  the  test.  The  maximum  voltage  23.76  mV  was  obtained  when 
applied  40  K  differential  temperatures.  We  studied  effect  of  geometric  design  by  using  a  TEG  with 
constant  cross  section  area  of  3x3  mm2,  but  variable  heights.  We  found  that,  when  increase  height,  output 
voltage  also  increase  in  linearly  fashion  until  reach  a  specific  height,  output  voltage  was  saturated.  After 
the  saturated  point,  stretched  more  height  do  not  gain  more  explicit  output  power.  This  information 
suggested  that  there  is  some  height  which  produced  optimal  output  voltage. 
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